#include <stdio.h>
#include <cstdlib>
/*#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
  */
     
int main (void){
    const gsl_rng_type * T;
    gsl_rng * r;
     
    int i, n = 10;
    double mu = 3.0;
     
    /* create a generator chosen by the 
    environment variable GSL_RNG_TYPE */
     
    gsl_rng_env_setup();
     
    T = gsl_rng_default;
    r = gsl_rng_alloc (T);
     
    /* print n random variates chosen from 
       the poisson distribution with mean 
        parameter mu */
     
    for (i = 0; i < n; i++) {
           unsigned int k = gsl_ran_poisson (r, mu);
           printf (" %u", k);
     }
     
       printf ("\n");
       gsl_rng_free (r);
       return 0;
     }
